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A WEIBULL  MODEL  TO  ESTIMATE  RESIDUAL  AND 
CRITICAL  VELOCITIES  FOR  TARGET  PENETRATIONS 


D.  CLARK 
L.  CROW 
J.  SPERRAZZA 


1.  INTRODUCTION 


A number  of  procedures  have  been  proposed  for  estimating  from  test 
data  the  functional  relationship  between  V , the  residual  velocity  of  a 

projectile  after  penetration  of  a target  and  V , its  striking  velocity. 

Examples  of  these  are  models  based  on  the  hyperbolic  [1]  and  exponential 
[4]  relationships 


/2  * AV2  + B and  V = V - V e8^1'  Vs/Vc\ 
r s r s c 


respectively.  In  these  models,  V , the  critical  velocity  is  defined  as 
the  Vs  intercept  when  V = 0. 

Another  penetration  prediction  model  is  the  Johnson  equation  [2] 

vr  = (Vs  - Vc)[eK4CVs  ‘ Vc)K5-i] 


where 


v = ^(sec  0-l)[eK2(A/M  x 103]K3_i] 


and  6 is  striking  obliquity  angle  in  degrees,  A is  fragment  presented 

area  in  cm2,  M is  fragment  weight  in  grams  and  Kj,  K2,  K3,  Kt,,  and  K5 

are  constants  to  be  estimated  from  the  data.  For  this  model  the  critical 
velocity  is  not  directly  estimated  from  test  data.  Instead,  the 
critical  velocity  is  determined  by  using  the  estimated  Vc's  from  an 

empirical  model,  such  as  the  hyperbolic,  fitted  to  several  sets  of 
data  to  estimate  the  constants  Kj,  K2  and  K3. 

For  meaningful  applications  of  these  prediction  models  in  ballistic 
studies,  it  is  important  that  they  re?.1 istically  represent  the  relation- 
ship between  striking  and  residual  velocity  and  provide  adequate  estimates 


of  the  critical  velocity  V . Various  applications  of  the  penetration 
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they  do  not  sufficiently  describe  this  relationship. 

In  this  report  we  propose  a procedure  based  on  the  versatile 
three-parameter  Weibull  distribution  function  for  estimating  the 
relationship  between  the  residual  and  striking  velocities  of  a projectile 
from  test  data.  This  model  has  many  shapes,  as  illustrated  in  the  next 
section,  which  should  make  it  useful  for  fitting  ballistic  data  over  a 
wide  range  of  firing  conditions  for  various  types  of  projectiles.  In 
Section  3 we  discuss  nonlinear  estimation  procedures  for  fitting  this 
model  to  test  data.  These  procedures  utilize  striking  velocities  with 
zero  residuals  to  help  estimate  the  three  unknown  parameters,  including 
the  critical  velocity.  In  Appendix  A we  list  a computer  program  for 
estimating  these  parameters  and  illustrate  its  use  by  fitting  the  Weibull 
model  to  several  sets  of  penetration  data.  In  Section  4 we  compare  the 
Weibull  and  hyperbolic  models  on  eight  sets  of  penetration  data. 

2.  THE  WEIBULL  MODEL 


Consider  a continuous  functional  relationship  = G(Vs),  between 

striking  velocity  and  residual  velocity,  with  other  impact  conditions 
fixed.  We  assume  that  the  function  G describing  this  relationship 
satisfies  the  following  conditions: 


(i) 

G(V  ) - 0 

0 < V 

— s 

(ii) 

G\Vs)  > 0 

V > V 

S i 

(iii) 

Vs  > GCVs) 

V > V 

S i 

(iv) 

Vs  - G(Vs)  - 0 

8 

t 

t/i 

> 

In  practice,  we  are  generally  concerned  with  determining  G over  a finite 
range  of  Vg  before  any  significant  fragmentation  of  the  projectile  occurs. 

However,  in  the  formulation  of  the  problem,  we  assume  the  asymptotic 
property  of  condition  (iv) . 


In  the  present  paper  we  model  the  functional  relationship  as 


V < V 
s — c 


G(Vs)  - 


(2.1) 


V [i-e_X^Vs"Vc^B]  V > V 

s L 1 sc 


where  A >0,  3 > 0 and  V >0.  Observe  that  G satisfies  conditions 

c 

(i)  - (iv)  and  has  (for  Vg  >_ V ) the  form  G(x)  = xF(x)  where 
A f i ® 

F(x)  =l-e  ^ ,x>n,  is  the  three-parameter  Weibull  distribution 

function  with  scale  parameter  A,  shape  parameter  3 and  location  parameter  n. 
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In  Figure  1 we  show  some  of  the  many  shapes  the  function  G(x)  can 
assume  for  various  values  of  X and  6 when  n a 100.  Condition  (iv) 
implies,  of  course,  that  the  line  y ■ x is  an  asymptote  of  G(x). 

However,  the  second  derivative  G"(x)  is  not  necessarily  less  than  zero 
for  all  x.  Hence,  as  shown  in  Figure  Id,  G(x)  may  actually  move  away 
from  the  line  y « x over  a finite  range  of  x,  for  certain  values  of 
X and  $.  However,  the  sign  of  G"(x)  will  eventually  change  and  G(x) 
will  approach  the  line  y * x asymptotically.  This  property  of  G(x) 
is  one  of  the  characteristics  which  makes  it  a versatile  model  for  fitting 
penetration  data. 

3.  ESTIMATION  PROCEDURES 

The  three  unknown  parameters  X,  $ and  in  the  model  given  by  (2,1) 

can  be  estimated  by  the  use  of  a nonlinear  programming  algorithm.  A 
program,  given  in  Appendix  A,  utilizing  Marquardt's  method  [5]  for 
nonlinear  least  squares  has  been  developed  by  AMSAA  for  this  application. 

Estimates  of  X,  0 and  V are  determined  as  values  X,  6 and  V , respectively, 
which  minimize  the  root  mean  square  error. 

ERMS  =^T[Vr  - GCVg)]2^1^, 

where  the  summation  is  taken  over  all  N pairs  (Vs>  V ) such  that 


The  estimates  X and  6 are  determined  by  this  procedure  with  the 
constraint  that  they  be  greater  than  zero.  Furthermore,  to  utilize 
the  information  associated  with  observed  residual  velocities  that  are 
zero,  we  perform  "his  optimization  with  the  additional  constraint  that 
a _<  V£  < b,  where  a and  b are  inputs  to  the  program. 

Marquardt's  algorithm,  which  is  an  unconstrained  optimization 
technique,  has  been  modified  to  accommodate  the  constraints  on  X,  6 
and  V . In  most  cases  the  algorithm  converges  to  the  solution  within 

very  few  steps.  If  satisfactory  results  are  not  obtained  by  use  of  this 
method,  one  may  elect  to  use  another  scheme  such  as  those  programmed  by 
Wortman  [6] . 

Regardless  of  the  method  employed  to  determine  optimum  values  of 
the  parameters,  preliminary  estimates  X°,  fj°,  and  V°  must  be  established. 

Discretion  should  be  exercised  in  order  to  determine  the  interval  [a,h) 

in  which  V is  constrained  to  lie.  It  is  often  feasible  to  choose  a to  be 
c 


11 


the  highest  striking  velocity  for  which  the  projectile  did  r.ot  perforate 
the  target  and  b chosen  to  equal  the  lowest  striking  velocity  for  which 
perforation  was  achieved.  Once  this  interval  has  been  determined  an 
initial  estimate  V“  may  be  arbitrarily  chosen  within  the  interval. 

With  this  choice  of  V“  the  equation  (2.1)  can  be  linearized  to  yield 

C 

Vs  \ 

-^-J  « ln(\)  + 0in(V  -V  ). 

s r J 

In  most  cases  initial  estimates  X°  and  e*  within  t.he  feasible  region  may 
be  then  obtained  by  using  linear  regression. 

In  Appendix  A we  illustrate  the  application  of  these  procedures 
and  the  versatility  of  the  model  using  several  sets  of  penetration  data. 

4.  COMPARISON  OF  WE I BULL  AND  HYPERBOLIC  MODELS 

In  a recent  study  [3]  published  by  the  USA  Ballistic  Research  Laboratories 

the  hyperbolic  model  ■ AV^  + B was  used  to  analyze  the  residual 

velocity  of  right  circular  cylinders  after  perforating  doron  body  armor 
material.  The  cylinders  considered  in  this  study  were  made  from  01  Tool 
Steel  and  heat  treated  to  a hardness  of  R£  29  + 2.  The  hyperbolic  model  was 

fitted  to  eight  sets  of  penetration  data  from  2,  4,  16  and  64  grain  cylinders 
at  0*  and  45°  obliquity. 

We  fitted  the  Weibull  model  to  the  eight  sets  of  data  considered 
in  the  above  report  and  compared  the  results  to  those  obtained  from 
the  hyperbolic  model.  These  comparisons,  with  the  fitted  curves,  are 
given  in  Tables  1-8  and  Figures  29,  respectively. 

In  four  cases  the  ERMS's  were  slightly  different  for  the  two  models 
(the  hyperbolic  ERMS's  being  lower  in  three  of  'ese  cases),  while  in 
the  other  four  cases  the  ERMS's  for  the  Weibull  ..iodel  were  significantly 
lower.  Also,  note  that  the  hyperbolic  model,  in  several  cases  estimated 
Vc  lower  than  what  one  would  expect  based  on  the  data.  For  example, 

consider  the  case  of  2 grain  steel  at  45  degrees  obliquity.  The  hyper- 
bolic model  estimated  V ac  530  m/s,  but  the  experimental  data  had  non- 
perforations for  striking  velocities  as  high  as  629  m/s  and  perforations 
for  striking  velocities  only  as  low  as  605  m/s.  The  Weibull  model  esti- 
mate of  V in  this  case  is  607  m/s, 
c 

5 . CONCLUSIONS 

The  Weibull  model  provides  a versatile  form  which  may  be  used  to 
represent  the  relationship  between  the  striking  velocity  and  t.he  resid- 
ual velocity  of  a projectile  fired  into  some  target  material.  This 
model  compares  favorably  with  others  proposed  for  this  application. 
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In  particular,  it  provides  more  satisfactory  estimates  of  the  critical 
velocity. 

< ' - • , 

To  date,  this  model  hps  been  used  only  to  fit  individual  sets  of 
data  for  which  all  conditions  other  than  the  striking  velocity  are 
held  constant.  No  attempt  has  been  made  to  physically  interpret  the 
values  of  the  parameter  estimates  or  to  interpolate  between  test  con- 
ditions. A model  with  those  capabilities  is  needed  for  use  in  vulner- 
ability and  effectiveness  models. 
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.0 

.0 

125.0 

44.0 

145.0 

199.0 

316.0 

376.0 

397.0 

443.0 

555.0 
5P1.0 

663.0 

736.0 

776.0 

833.0 

894.0 

974.0 

1002.0 

1031.0 


STRIKING 

VELOCITY 

M/S 

308.0 

423.0 

431.0 

446.0 

424.0 

451.0 

466.0 

486.0 

537.0 

597.0 

622.0 

662.0 

773.0 

807.0 

879.0 

975.0 

1027.0 

1067.0 

1143.0 

1247.0 

1290.0 

1306.0 


RESIDUAL 

VELOCITY 

M/S 

.0 

.0 

.0 

.0 

75.0 

87.0 

141.0 

110.0 

242.0 

320.0 

360.0 

415.0 

535.0 

569.0 

638.0 

721.0 

768.0 

794.0 

851.0 

945.0 

969.0 

990.0 
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TABLE  3 


16  CRAIN  STEEL  RT  CIRCULAR  CYLINDER  INTO 
STANDARD  OORON  AT  0 OEGREES  OBLIQUITY 


2 2 

HYPERBOLIC  MODEL  VR  * -58815.3  ♦ .7782  VS 

CRITICAL  VELOCITY  = 274.0  M/S 

ERROR-RMS  * 15.3  M/S 


.313 

HEIBULL  MODEL  VR/VS  = l-EXPI-  .2440IVS-  300.1) 

CRITICAL  VELOCITY  = 309.1  M/S 

ERROR-RMS  * 16.1  M/S 


EXPERIMENTAL  DATA 


STRIKING 

RESIDUAL 

VELOCITY 

VELOCITY 

M/S 

M/S 

282.0 

.0 

287.0 

.0 

322.0 

120.0 

397.0 

236.0 

537.0 

404.0 

693.0 

562.0 

874.0 

747.0 

953.0 

808.0 

1109.0 

940.0 

1251.0 

1069.0 

STRIKING 

RESIDUAL 

VELOCITY 

VELOCITY 

M/S 

M/S 

284.0 

.0 

312.0 

96.0 

349.0 

170.0 

477.0 

333.0 

615.0 

482.0 

773.0 

656.0 

881.0 

756.0 

1037.0 

882.0 

1195.0 

1024.0 
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TABLE  4 


64  GRAIN  STEEL  RT  CIRCULAR  CYLINDER  INTO 
STANDARD  OORON  AT  0 DEGREES  OBLIQUITY 


2 2 
HYPERBOLIC  MODEL  VR  * -31554.4  ♦ .8533  VS 

CRITICAL  VELOCITY  = 192.3  M/S 

ERROR-RMS  * IT. 6 M/S 


.2623 2 

NEI BULL  MODEL  VR/VS  = 1-EXPI-  .4322IVS-  245.0) 

CRITICAL  VELOCITY  * 245.0  M/S 

ERROR-RMS  * 16. 6 M/S 


EXPERIMENTAL  CATA 


STRIKING 

VELOCITY 

M/S 

193.0 

223.0 

249.0 

277.0 

278.0 

308.0 

314.0 

384.0 

551.0 

633.0 

635.0 

809.0 

989.0 


RESIDUAL 

VELOCITY 

M/S 

.0 

.0 

121.0 

163.0 

160.0 

240.0 

260.0 

293.0 

479.0 

558.0 

560.0 

728.0 

899.0 


STRIKING 

RESIDUAL 

VELOCITY 

VELOCITY 

M/S 

M/S 

208.0 

.0 

228.0 

.0 

250.0 

118.0 

277.0 

200.0 

300.0 

200.0 

309.0 

203.0 

326.0 

222.0 

461.0 

418.0 

554.0 

488.0 

635.0 

543.0 

707.0 

638.0 

927.0 

843.0 

1153.0 

1043.0 
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TABLE  5 


2 GRAIN  STEEL  RT  CIRCULAR  CYLINDER  INTO 
STANDARD  OORON  AT  45  DEGREES  OBLIQUITY 


2 2 
HYPERBOLIC  MOOEL  VR  *-114876.9  ♦ .4091  VS 

CRITICAL  VELOCITY  = 529.9  M/S 

ERROR-RMS  = 51.3  M/S 


.25822 

MEIBULL  MODEL  VR/VS  = 1-EXP(-  .1656(VS-  606.91  ) 

CRITICAL  VELOCITY  = 606.9  M/S 

ERROR-RMS  = 40.7  m/s 


EXPERIMENTAL 

DATA 

STRIKING 

RESIDUAL 

STRIKING 

RESIDUAL 

VELOCITY 

VELOCITY 

VELOCITY 

VELOCITY 

H/S 

M/S 

M/S 

M/S 

544.0 

.0 

561.0 

.0 

577.0 

.0 

587.0 

.0 

587.0 

.0 

600.0 

.0 

601.0 

.0 

609.0 

.0 

615.0 

.0 

617.0 

.0 

621.0 

.0 

629.0 

.0 

605.0 

72.0 

607.0 

51.0 

620.0 

140.0 

628.0 

170.0 

632.0 

194.0 

635.0 

235.0 

636.0 

183.0 

657.0 

192.0 

671.0 

208.0 

678.0 

184.0 

680.0 

217.0 

714.0 

361.0 

740.0 

378.0 

781.0 

418.  0 

789.0 

4r6 . 0 

798.0 

417.0 

855.0 

413.0 

863.0 

529.0 

877.0 

421.0 

901.  C 

525.0 

931.0 

497,0 

960.0 

526.0 

1018.0 

574.0 

1067.0 

619.0 

1078.0 

595.0 

1172.0 

644.0 

1209.0 

696.0 

1221.0 

679.0 

1235.0 

698.0 

1341.0 

780.0 

1414.0 

839.0 

1457.0 

867.0 
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TABLE  6 


4 GRAIN  STEEL  RT  CIRCULAR  CYLINDER  INTO 
STANDARD  OORON  AT  45  OEGREES  OBLIQUITY 


? 2 
HYPERBOLIC  MODEL  VR  = -72019.2  ♦ .5039  VS 

CRITICAL  VELOCITY  = 378.1  M/S 

ERROR-RMS  * 55.2  M/S 


.20376 

MEI BULL  MODEL  VR/VS  = l-EXP(-  .3015IVS-  494.3)  ) 

CRITICAL  VELOCITY  = 494.3  M/5 

ERROR-RMS  * 40.9  M/S 


EXPERIMENTAL 

DATA 

STRIKING 

RESIDUAL 

STRIKING 

RESIDUAL 

VELOCITY 

VELOCITY 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

423.0 

.0 

437.0 

.0 

454.0 

.0 

459.0 

.0 

467.0 

.0 

474.0 

.0 

475.0 

.0 

480.0 

.0 

481.0 

.0 

483.0 

.0 

483.0 

96.0 

495.0 

134.0 

497.0 

123.0 

505.0 

148.0 

536.0 

204.0 

575,0 

275.0 

629.0 

366.0 

671.0 

411.0 

688.0 

428.0 

717.0 

459.0 

799.0 

544.0 

840.0 

585.0 

919.0 

646.0 

1002.0 

694.0 

1171.0 

787.0 

1285.0 

866.0 

1299.0 

872.0 

1348.0 

893.0 

1348.0 

912.0 
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STRIKING  VELOCITY  M/S 

IRCULRR  CYLINDER  INTO  STRNDRRD  DORON  RT  45  DEGREES  OBLIQUITY 

FIGURE  7 


TABLE  7 


16  GRAIN  STEEL  RT  CIRCULAR  CYLINDER  INTO 
STANOARO  OORON  AT  45  DEGREES  OBLIQUITY 


? 2 t 

HYPERBOLIC  MODEL  VR  = -50639.3  ♦ ,6889  VS  j 

CRITICAL  VELOCITY  = 271.1  M/S  | 

ERROR-RMS  = 35.9  M/S  i 

i 

i 


.28592 


BULL  MODEL 

VR/VS  = 1-EX°( 

- .25151 VS- 

323.31 

CRITICAL  VELOC 

ITY  = 323. 

3 m/s 

ERROR-RMS  - 

26. 

6 M/S 

EXPERIMENTAL 

DATA 

STRIKING 

RESIDUAL 

STRIK INC 

RESIDUAL 

VELOCITY 

VELOCITY 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

302.0 

.0 

304.0 

.0 

304.0 

.0 

313.0 

.0 

324.0 

.0 

315.0 

38.0 

324.0 

78.0 

327.0 

90.0 

335.0 

96.0 

335.0 

00.  C 

348.0 

111.0 

370.0 

169.0 

408.0 

234.0 

433.0 

278.0 

486.0 

367.0 

508.0 

361.0 

545.0 

393.0 

567.0 

421.0 

587.0 

441  .0 

665.0 

518.0 

699.0 

549.0 

79  7.0 

639.0 

859.0 

693.0 

905,0 

734.0 

955.0 

772.0 

964.0 

778.0 

979.0 

791.0 

1005.0 

811.0 

1026.0 

834.0 

1138.0 

919.0 

1159.0 

935.0 

1207.0 

944.0 

1227.0 

984.0 

1234.0 

990.0 

j' 

i 
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TABLE  8 

64  GRAIN  STEEL  RT  CIRCULAR  CYLINDER  INTO 
STANOARD  OCRON  AT  45  DEGREES  OBLIQUITY 


2 2 
HYPERBOLIC  MODEL  VR  « -15795.6  ♦ .7950  VS 

CRITICAL  VELOCITY  = 141.0  M/S 

ERROR-RMS  = 38.2  */S 


.18142 


BULL  MODEL 

VR/VS  = l-pXP( 

- .66281  VS- 

240.0) 

CRITICAL  VELOCITY  = 240. 

0 M/S 

FRROR-RMS  = 

31. 

8 M/S 

EXPERIMENTAL 

DATA 

STRIKING 

RESIDUAL 

STRIKING 

RESIDUAL 

VELOCITY 

VELOCITY 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

211.0 

.0 

212.0 

.0 

213.0 

.0 

215.0 

.0 

225.0 

77.0 

233.0 

64.0 

243.0 

120.0 

244.0 

116.0 

271.0 

162.0 

275.0 

165.0 

334.0 

290.0 

379.0 

290.0 

383.0 

362.0 

451.0 

412.0 

535.0 

496.0 

626.0 

576.0 

674.0 

588.0 

733.0 

644.0 

790.0 

671.0 

792.0 

697.0 

852.0 

755.0 

904.0 

803,0 

967.0 

859.0 

1013.0 

902.0 

1067.0 

948.0 

1114.0 

987.0 

1159.0 

1028. 0 

1203.0 

1042.0 
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STRIKING  VELOCITY  M/S 

6H  GRAIN  STEEL  RT  CIRCULAR  CYLINDER  INTO  STANDARD  DORON  AT  45  DEGREES  OBLIQUITY 
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Appendix  A 


A FORTRAN  computer  program  has  been  prepared  to  provide  estimates 
of  the  parameters  for  the  Weibull  model  described  in  this  report.  The 
program  was  written  for  use  on  the  Ballistic  Research  Laboratories 
computers,  BRLESC  I and  II.  The  program  utilizes  16K  memory  locations. 
The  CALCOMP  plotter  package  is  used  to  provide  plots  of  the  data  and 
fitted  curves. 

The  required  input  is  as  follows: 


Columns 

Card  1 

1-79 

Alphanumeric  string  used  for  title 

80 

The  symbol 

Card  2 

1-10 

Minimum  value  for  critical  velocity 

11-20 

Maximum  value  for  critical  velocity 

Card  3 

1-10 

Striking  velocity 

11-20 

Residual  velocity 

Repeat  Card  3 as  necessary  to  complete  the  data  set.  Follow  the 
last  Card  3 with  a blank  card. 

After  execution  of  the  computation  for  a data  set  the  program 
returns  to  read  card  1 for  another  set. 

A listing  of  the  source  program  is  provided  on  the  following  pages. 

To  provide  several  samples  of  output  and  also  to  demonstrate  the  flex- 
ibility of  the  Weibull  model  the  output  from  several  data  sets  is  also 
included.  These  data  were  obtained  from  firings  of  a 90  percent  tungsten 
spheriod  weighing  approximately  7 grains.  The  firing  conditions  include 
two  target  materials,  two  thicknesses,  and  two  angles  of  obliquity. 


PRCGRAW  WBLRV 

CCPMCN  /H PX/  HHAX,  HWIN 

DIMENSION  VS(ICO),  VR(IOO),  ZU00«2t*  A(2,3),  C(3),  RtlOOl, 
1 AF ( 100 1 f SIG<3),  T(3),  m(8|,  F(IOO),  TXOIf  TV  131 
OIFENSION  XZ (20) , YZ (20) 

EXTERNAL  VREF 

DATA  T * ( Lit  TX(2),  YY<11,  TY(2)  /10HSTR1KING  V,  8HEL0CITV  , 
1IOLAL  , 9HVEL0CI TV  / , TX(3I,  TY(3)  /4HH/S>,  4HfVS>/ 

IAN*  1 

l READ  (5,111  (ITUI), 1*1,8) 

WRITE  (6,13)  (I?l(I),I«l,8l 
READ  (5,16)  HMIN.HItAX 
WRITE  (6,18)  HP!  N,Hf*AX 
F = 1 
NZ*C 

l READ  (5,12)  V$(P),VR(M) 

IF  ( VS(M).lE.O, ) GO  TO  4 
IF  (VR(P).LE.O. ) GO  TO  3 
IF  < VS(P).LE.HKIN)  GO  TO  3 
P-P+l 
GC  TO  2 
) N Z = N Z ♦ l 

XZ(NZI*VS(P) 

YZ(NZ  )*VR  IP) 

GO  TC  2 
F*P-1 

IF  (NZ.GT.C.)  CALL  SORTXY  (XZ,YZ,NZt 
CALL  SCRTXY  (VS.VR.H) 

(-‘*FM  IN 
DC  5 K*I  ,P 

F(K)*ALOG(ALCG(VS(K)/(VS(R)-VR(K)) ) ) 

Z(K,1)*1. 

Z(K.2I*ALCG(VS(K)-H) 

5 CONTINUE 

IF  (M.GT.2)  GO  TO  6 
WRITE  (6,17)  ft 
GO  TO  1 

6 CALL  GENL SQ  ( Z , 100 ,F ,M, A ,2 ,2 , C ,R , AF ,FRMS , S IG,T ,DET, 1 ) 

DO  7 1=1, P 

F ( I ) « VR  ( I ) 

Z(I,1)*VS(I) 

7 Z(I,2)-VR(I) 

C ( 3 ) *C  ( 2 ) 

C ( 2 ) *EXP ( C ( 1 1 ) 

C<  1)  =*H 

WRITE  (6,15)  (C(I), 1-1,3) 

CALL  PRQLS  ( VREP.F ,Z ,C ,2 ,P, 100 , 3, ERMS, S IG,R, AF, , IE-5 , . IE-3 ) 
8-C ( 3 ) 

XL  *C ( 2 ) 

H*C(  1) 

WRITE  (6,20  H,XL,B,ERMS 
WRITE  (6,21) 

IF  (NZ.GT.O)  WRITE  (6,15)  ( XZ ( I ) , YZ (!) , I* l,NZ ) 

DO  8 1*1, P 

8 WRITE  (6,1A)  VS(II,VR(I),AF(It,R(!) 

VPAX*  15CC. 

YS-YMAX/6. 

XOR-O. 

XHAX-X0R+8.PYS 

CALL  PLTVG  ( VS, VR ,M»XOR ,0. ,V$ , YS,TX,TY, 2,5, INN, ITL  ) 

IF  (NZ.GT.O)  CALL  PITCCO  (2 ,5 ,XZ ( 1 ) „YZ ( 1 ) ,NZ ) 

VS ( 1 ) *H 
VR( l)-0. 

34 


H, XL, B, ERRS 
WRITE  (6,151 


(XZ(I),YZ(I),I*l,NZ ) 


VS(I1,VR(I),AF(I),R(!) 


WBLRV  l 
WBLRV  2 
WBLRV  3 
WBLRV  4 
WBLRV  5 
WBLRV  6 
10HRESWBLRV  70 
WBLRV  80 
WBLRV  5 
WBLRVIOR 
WBLRV 1 1 W 
WBLRV12R 
WBLRV13W 
WBLRV14 
WBLRV15 
WBLRV16R 
WBLRV17 
WBLRV1B 
WBLRV19 
WBLRV20 
WBLRV21 
WBLRV22 
WBLRV23 
WBLRV24 
WBLRV25 
WBLRV26 
WBLRV27 
WBLRV28 
WBLRV29 
WBLRV30 
WBLRV31 
WBLRV32 
WBLRV33 
WBLRV34 
WBLRV35 
WBLRV  36W 
WBLRV37 
WBLRV38 
WBLRV39 
WBLRV40 
WBLRV41 
WBLRV42 
WBLRV43 
WBLRV44 
WBLRV45 
WBLRV46W 
W8LRV47 
WBLRV48 
WBLRV49 
WBLRV50 
W0LRV51W 
WBLRV52W 
WBLRV53W 
WBLRV54 
WRLRV55H 
WBLRV56 
WBLRV57 
WBLRV58 
WBLRV59 
WBLRV60 
WBLRV61 
WBLRV62 
WBLRV63 


( 


1C 


DX«.C25*YS 

1-1 

VS»V«VRUUOX 

vsi-vsm+ox 

IF  ( V ST . GT . X K A X ) GO  TO  10 
I » I ♦ 1 

vsi  u«vst 

VRI  I )»VSII)MI.-EXP<-XL*IVS<1  )-H)**B) ) 

VSAV«VR< l t 

IF  lt.GT.40)  0X».25*YS 
IF  IVSAV.LT.VMAX)  GO  TO  9 

CALI  PLTCCO  (l,,V$m,VRm,I.O,XOR,XMAX.O.,YRAX) 

INNSC 

GO  TC  1 


11  FORMAT  ( P A 1 C ) 

12  FORMAT  ( 2F1C.0) 

13  FORMAT  I 1H1.30X.4AIO/33X.3A10.A9/) 

14  FORMAT  ll9X,4F14.1t 

15  FORMAT  ( 1SX.2F14.L) 

It  FORMAT  12F1C.5) 

17  FORMAT  15HCCMY.I3.25H  DATA  PC1NTS--FIT  OMITTED) 
ie  FORMAT  I 2CX , 53MCR I T I CAL  VELOCITY  IS  CONSTRAINED  TO  THE 
1CM,F7.1/2CX, 3HTC  .F7.1.5H  M/S.//) 

19  FORMAT  I 2CX  . 1 3HW6 I BULL  MODEL » flX , 6HVC  M/S > 6X , 24HL AMBC A 
l ERMS//2CX, 17HINI TIAL  ESTIMATE S ,F 10 . 1 , 3X , 2F 10 .6/ ) 

2C  FORMAT  I 2CX.15HFINAL  E ST  I MATES ,F 12 . 1 . 3X , 2F 10 .6, F7 . 1// / ) 
21  FORMAT  ( 2tX, A 9H STRIKING  RESIDUAL  APPROXIMATION 

16X  » 22HVELCC  I T Y VELOCITY/1  P Xt  A < 1 1 X • THi'^S  ) / ) 

END 

SUBROUTINE  VREM  M> , X ,F  , P , IC , M , K ) 

COMMON  /HMX/  HMAXfHMIN 
DIMENSION  BIK)  , XtM) . P(K) 

IF  (B(l).GT.HMAX)  Bill *HM  A X 
IF  (Bllt.LT.HMIN)  BU)«HMIK 

xp»x<  it-pui 

IF  I XP.GT *0. ) GO  TO  l 
f-C. 

GO  TC  2 

1 XPL=ALOG(XP) 

XPP  = EXP(M3)*XPL) 

EX“EXP(”fH2)*XPB) 

F*X< 1>*( t.-EX) 

2 IF(IC.NE.C)  RETURN 
IFIXP.GT.OGC  TO  3 
P<1I«0. 

P ( 2 ) *C. 

PI 3>«C. 

RETURN 

3 Pill— EX*XU)*B(2)*B(3)*XPB/XP 
PI2)«X<1)*XPB*EX 
P(3)*PI2)*BI2)*XPl 

RETURN 

ENC 

SUBROUTINE  MRQLS  I FORM , Y , X ,R ,M , N, NM AX , K , FRMS , SE , R , Ff T AU 
D1MFNSI0N  YIN).  XINMAX.M),  BIK),  SEIK),  R I N ) » FIN) 
DIMENSION  A 1 20. 20 ) . SAI20.20),  PI20),  VI10).  GI20),  SGI 
GNL* 1C. 

IC  T “0 
Xl«.Cl 


WBLRVtA 
WBLRV65 
WBLRV66 
WBLRV67 
WBLRV68 
WBLRV69 
WBLRV70 
WBLRV71 
WBLRV72 
WBLRV73 
WBLRV7A 
MBLRV75 
WBLRV76 
HBLRV77 
MBLRV78 
WBLRV79 
WBLRV80 
WBLRV61 
MBLRV62 
MBLRV83 
MBLRVUA 
HBLRV85 
INTERVAL  FRWBLR V86 
WBLRV87 
BETA  WBLRV88 
MBLRV89 
MBLRV90 
FRR0R/2WBLRV91 


.EPS) 

20) 


UBLRV92 

MBLRV93- 

VREM 

1 

VREM 

2 

VREM 

3 

VREM 

A 

VREM 

5 

VREM 

6 

VREM 

7 

VREM 

8 

VREM 

9 

VREM 

10 

VREM 

11 

VREM 

12 

VREM 

13 

VREM 

1* 

VREM 

15 

VREM 

16 

VREM 

17 

VREM 

18 

VREM 

19 

VREM 

20 

VREM 

21 

VREM 

22 

VREM 

23 

VREM 

24 

MRQLS 

1 

MRQLS 

2 

MRQLS 

3 

MRQLS 

i A 

MRQLS 

5 

MRQLS 

6 

1 I c * c 

STEP*!* 

ICT* ICT* l 
PH  I *0  « 

00  2 1*1.20 
G< I t-G. 

0C  2 J*1 .20 

2 AU»J)*0. 

DO  5 1*1. N 
OC  3 J*l.H 

3 V(JJ*XU.J) 

CALL  FORM  (P.V.Fm.P.K.P.KI 

mii>yiiKin 

PHI-PHI+RI! >**2 

CC  A J-l.K 

GUI>GUl4R(tl*PIJI 

DO  A l*j,;< 

A A(J,L)*A1J,L)+P(J) *P (LI 

5 CONTINUE 

ERM$*SQRT (PHI /FLOAT  IN) > 

OC  6 1*1, K 
SEm*SQRTUU,!»» 

6 G<  1)*G(I  l/SEUI 
OC  7 1*1, X 

OC  7 J*I »K 

AU,jt*AU,j)/(SEm*SEun 
IF  1J.GT.I1  A(J.I1*AII,J1 

7 CONTINUE 
00  8 1-1. K 
SG{!)«G(I1 
DO  8 J-l.K 

8 SA(I.J)-Afl.J) 

1C  * 1 

XL  P*  XL/GNU 

9 DC  1C  1*1, K 

10  A1 1 , I )-A ( I , l 1 + XLH 

CALL  HAT  I NV  1 A ,K ,G ,20 , l ,DET 1 
IF  ( DET . "C  »0.  ) GO  TO  17 

11  00  12  I-l.K 

12  GUS-eil  )4STEP*IG<ll/SEim 
PH  IL-O* 

OC  1A  I-l.N 
00  13  J-l.H 

13  V1J1-XU.J) 

CALL  FORM  (G.V.Ftll.P.IC.McM 
RUl*vm-F(II 
1 A PHIL*PH1L*R(I )♦*£ 

GO  TO  115,10).  IC 

15  IF  1PHIL-PHI)  16,16,17 

16  XL-XLH 
GO  TO  25 

17  00  18  1*1, K 
Gdl-SGd) 

00  18  J«l,K 

ie  A 1 1,  J ) *SA  C I , J) 

IC-2 
XLH*XL 
GO  TC  3 

19  IF  (PH1L-PHI ) 25,25,20 

20  IF  1STEP.LT. 1.)  GO  TO  ?A 


MRQtS  7 
MRQtS  8 
MRQLS  9 
MRQL310 
MRQtSll 
MRQLS  12 
MRQLS 13 

MRQLS 1A 

MRQLS 15 

MRQLS 16 

MRQLS 17 

MRQLS18 

MRQLS 19 

MRQLS20 

MRQLS21 

MRQLS22 

MRQLS23 

MRQLS2A 

MRQLS25 

MRQLS26 

MRQLS27 

MRQLS28 

MRQLS 29 

MRQLS30 

MRQLS31 

MRQLS32 

MRQLS33 

MRQLS3A 

MRQLS35 

MRQLS36 

MRQLS 37 

MRQLS 38 

MRQLS39 

MRQLSAO 

NRQLSA1 

MRQLSA2 

MRQLSA3 

MRQLS AA 

MRQLSA5 

MRQLSA6 

MRQLSA7 

MRQLSA8 

MRQLSA9 

MRQLS50 

MRQLS51 

MRQLS52 

MRQLS53 

MRQLS5A 

MRQLS55 

MRQLS56 

MRQLS57 

MRQLS58 

MRQIS59 

MRQLS60 

MRQLS61 

MRQLS62 

MRQLS63 

MR0LS6A 

MRQLS65 

MRQLS66 
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IF  ( XL.Gf  .GNL*GNU)  GO  TO  22 

21  XL*XL*ONL 
GC  1C  17 

22  GNCRH2*G. 

CNCBM2*C. 

COT*C. 

OC  23  1 * 1 
GI«SGU)*$E(J) 

Ol*G<  Il-BU  ) 

OOT«DCT»GI*OI 

GNCRH2*0\CRK?+GI*GI 

23  DNCRF2*r> CRP2*Cl*OI 
OCT*DOT/SC«T (GNCRP2*CNORM2) 

24  IF  (DCT.LE..7C71)  GC  TO  21 
STEP=.5*STEP 

GC  T 2 11 

25  CC  26  I - 1 » K 

IF  (ABMGU  )-BU)  )/<TAlJ«  ABSIRU  m.GT.EPSI  GO  TO  29 

26  CCMINUE 

27  ERFS=SQRT (PHIL/FLOAT ( N ) 1 
CC  28  1*1, K 

28  B I 1 ) *G ( 1 I 
RETURN 

29  DC  2G  1*1, K 
2C  fl(  II»G(n 

IF  nCT.LT. 50  CO  TO  1 
GO  TO  27 

31  FCRVAT  UFC,  1BHFAILEC  TO  C0KVFPG1) 

EMC 

SlIBR  COT  I NE  PLTVG  ( X , V , N , XCR  , Y0*  , XS  » YS  ,TX  , T V , P , NS  , INN,  I TL  ) 
CIPrN$lCN  X(N),  YIN),  TX<?),  TY<2),  815000),  XCC2),  YC(2), 
IF  I INN.NE.l ) GO  TO  1 
CALL  PLTCCR  (12. ,1,6(1)  ,R(5COO) 

IC  T = 2 

1 I * I C.  T ♦ 1 

ICT  *HCO( 1,3) 

IF  ncr.EC.C)  CALL  PLTCCP 
XB  = 1 . 

YB=I .49. 5AFLCAT I ICT) 

CALL  PLTCCS  < xe,YB,0. ,0. ,1. ,1. ) 

CC  1 1*1,2 

XP * 1 C • 5*FL0 A T ( I-l ) 

OC  2 J*l,2 
YP*P.*FLCAT(J-1) 

XC(1 )-XP 
XCC2)*XP 
YCI 1 )=YP-#25 
YC  ( 2 ) =Y<*4 . 2 5 

CALL  PLTCCO  u,o,xcm  ,VC<i)  ,2) 

XCI 1 )=<P-.25 
XC  ( 2 ) *XP-»  ,25 
YC ( 1 ) *YP 
YC(2)*VP 

2 CALL  PLTCCO  ( 1 , 0 , XC ( l ) , YC ( 1 ) , 2 ) 

3 COMINU  = 

UlS.TXU  ) ,3.  ,1.  ,4., .75) 

(.15»TY(1) fl.,0.,.9,3.) 
I.lOflTL(l),0.,l.,l.5,.B) 


CALL  PLTCCT 
CALL  PUCCT 
CALL  PLTCCT 
X8-XB+1.5 
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MRQIS67 
HRQLS68 
HRQLS69 
HRQLS7C 
MRQLS71 
HRQLS72 
MRQLS73 
HRQLS74 
KRQLS75 
MRQLS76 
HRQLS77 
MRQLS78 
MRQLS79 
MRQLS80 
MRQLSti 1 
MRQLS82 
MRQLS83 
MRQLS84 
HRQLS85 
MRQLS86 
MRQLS87 
MROLsea 
MRQLS89 
MR0LS90 
MRQLSS1 
MRQLS92 
MRQLS94 
MRQLS95 
HRQLS96 
MRQLS97- 
PLTVG  1 
ITLI8)  PLTVG  2 
PLTVG  3 
PLTVG  A 
PLTVG  5 
PLTVG  6 
PLTVG  7 
PLTVG  8 
PLTVG  9 
PLTVG  10 
PLTVG 1 1 
PLTVG  12 
PLTVG  13 
PLTVG  14 
PLTVG 15 
PLTVG16 
PLTVG17 
PLTVG  18 
PLTVG19 
PLTVG2C 
PLTVG21 
PLTVG22 
PLTVG23 
PLTVG24 
PLTVG25 
PLTVG26 
PLTVG27 
PLTVG28 
PLTVG29 
PLTVG30 


YB»YB*1.5 

PLTYG3I 

CALL  PLTCCS  ( XB»YB  »XOR ,YOR*X$  #YS > 

PLTYG32 

XMAX-XOR+B.*XS 

PLTYG33 

YMAX«Y0R46.*YS 

PLTYG34 

CALL  PLTCCA  ( XS . YS ,X0R ,XMAX , YOR * YMAX.O ) 

PLTWG35 

CALL  PLTCCO  < PfNS »X ( U »Y( I ) ,N ) 

PLTVG36 

CALL  LABELA  IXS,YS,XOR»XHAX .YOR.YMAX, 1. . I. ) 

PLTVG37 

RETURN 

PLTVG38 

ENC 

PLTYG39- 

♦ 

COMPILE  DISC, 6ENLSQ, ALL 

♦ 

COMPILE  CISCtSORTXYtALL 

1* 

COMPILE  DISC. LABELA. ALL 

c 

ENC 

\ 


i 


! 

i 


t. 


i 


« 
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TABLE  A-l 


.45  GM  TUNGSTEN  SPHEROID  INTO  .15  CM 
ALUMINUM  AT  0 DEGREES  OBLIQUITY 

CRITICAL  VELOCITY  IS  CONSTRAINED  TO  THE  INTERVAL  FROM  243.8 
TO  310.0  M/S. 


ME  I BULL  MODEL  VC  M/S  LAMBDA  BETA  ERMS 

INITIAL  ESTIMATES  243.8  ,246400  .328262 


ESTIMATES 

298.0 

.504504  .222161 

13.0 

STRIKING 

REST  DUAL 

approximation 

ERROR 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

243.8 

.0 

302.4 

153.0 

152.1 

1.0 

355.7 

246.3 

253.0 

-6.7 

363.3 

271.0 

262.0 

9.0 

373.4 

253.3 

273.4 

-20.1 

374.0 

766.4 

275.1 

-8.7 

399.9 

321.3 

30  2.2 

19.1 

422.8 

318.7 

325.9 

-7.  7 

976.9 

880.9 

862.8 

18.0 

1026.3 

905.9 

910.4 

-4.6 

1033.3 

921.7 

917.2 

4.5 

1036.6 

940.? 

9 20.4 

19.9 

1494.1 

1371.6 

1363.3 

8.3 

1501.4 

1353.0 

1370.4 

-17.4 

1547.2 

1400.6 

1414.9 

-14.3 

! 
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STRIKING  VE 

-H5  GM  TUNGSTEN  SPHEROID  INTO  .15  CM  fil 

FIGURE  A-l 


TABLE  A- 2 


.45  GM  TUNGSTEN  SPHEROID  INTO  .32  CM 
ALUMINUM  AT  0 DEGREES  OBLIQUITY 


CRITICAL  VELOCITY  IS  CONSTRAINED  TO  THF  INTERVAL  FROM  343.5 
TO  460.0  M/S, 


MEIBULL  MOOEL  VC  M/S 

INITIAL  ESTIMATES  343.5 


LAMBDA 

BETA 

ERMS 

.426603 

.214394 

ESTIMATES 

343.5 

.333039  .253525 

40 

STRIKING 

residual 

APPROXIMATION 

ERROR 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

343.5 

.0 

456.9 

379.5 

305.6 

73.9 

526.4 

307.8 

375.2 

~67.3 

623.3 

423.1 

468.0 

-44.9 

899.2 

718.7 

727.1 

-8.3 

915.3 

766.0 

742.2 

23.7 

920.8 

770.5 

747.4 

23.2 

1524.0 

1326.2 

1318.0 

8.1 

1531.6 

1322.2 

1325.3 

-3.  1 
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STRIKING  VELOCITY  M/S 

H5  GH  TUNGSTEN  SPHEROIO  INTO  -32  CM  RLUMINUM  RT  0 DEGREES  OBLTQUtTY 

FIGURE  A-2 


TABLE  A-2 


.45  GM  TUNGSTEN  SPHEROID  INTO  .32  CM 
ALUMINUM  AT  0 OEGREES  OBLIQUITY 

CRITICAL  VELOCITY  IS  CONSTRAINED  TO  THE  INTERVAL  FROM  343.5 
TO  460.0  M/S. 


ME  I BULL  MODEL 


VC  M/S 


LAMBDA  BETA  ERMS 


INITIAL  ESTIMATES  343.5 


.426603  .214394 


FINAL 


ESTIMATES 

343.5 

STRIKING 

RESIDUAL 

VELOCITY 

VELOCITY 

M/S 

M/S 

343.5 

.0 

456.9 

379.5 

526.4 

'07.8 

623.3 

423.1 

899. 2 

718.7 

915.3 

766.0 

920.8 

770.5 

1524.0 

1326.2 

1531.6 

1322.2 

.333039  .253525  40.7 


APPROXIMATION  ERROR 

M/S  M/S 

305.6  73.9 

375.2  -67.3 

46B.0  -44.9 

727.1  -8.3 

742.2  23.7 

747.4  23.2 

1318.0  8.1 

1325.3  -3.1 
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•FABLE  A-3 


.45  GM  TUNGSTEN  SPHEROID  INTO  .15  CM 
MILO  STEEL  AT  0 DEGREES  OBLIQUITY 

CRITICAL  VELOCITY  IS  CONSTRAINED  TO  THE  I NT  E I /AL  FROM  343 
TO  380.0  M/S. 


WEI  BULL  MODEL  VC  M/S 
INITIAL  ESTIMATES  343.5 
FINAL  ESTIMATES  358.1 


LAMBDA  BETA  ERMS 

.213101  .292740 

,291155  .244159  56.8 


STRIKING 

RESIDUAL 

VELOCITY 

Vc!.OCtTY 

M/S 

M/S 

281.6 

.0 

305.4 

.0 

318.5 

.0 

334.1 

.0 

343.5 

.0 

365.5 

138.4 

388.3 

183. R 

391.4 

204.2 

410.3 

208.2 

559.3 

370.9 

883.6 

711.4 

920.8 

710.2 

920.8 

696.8 

961.6 

780.0 

989.7 

755.3 

1026.3 

596. S 

1524.0 

1228.0 

1531.6 

1257.9 

APPROXIMATION  ERROR 

M/S  M/S 


138.0 

.3 

189.5 

-5.7 

194.1 

10.1 

219.3 

-11. 1 

366.1 

4.8 

653.2 

58.2 

686.0 

24.2 

686.0 

10.8 

722.2 

57.8 

747.0 

8.3 

779.5 

-183.0 

1226.2 

1.9 

1233.1 

24.8 
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TABLE  A-4 


.45  GM  TUNGSTEN  SPHEROID  INTO  .32  CM 
MILD  STEEL  AT  0 DEGREES  OBLIQUITY 


CRITICAL  VELOCITY 
TO  530.0  M/S. 

IS  CONSTRAINED 

TO  THE 

INTERVAL  FROM 

515 

WEI  BULL  MODEL 

VC  M/S 

LAMBDA 

BETA 

ERMS 

INITIAL  ESTIMATES 

515.7 

.108057 

.324210 

FINAL  ESTIMATES 

515.7 

•08949? 

.356297 

29.0 

STRIKING 

RESIDUAL 

APPROXIMATION 

ERROR 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

515.7 

.0 

520.9 

92.4 

77.5 

14.9 

698.9 

229.2 

210.3 

18.9 

603.5 

189.0 

215.1 

-26.2 

607.2 

239.0 

219.0 

20.0 

873.3 

389. 6 

451.2 

-61.6 

909.8 

511.1 

481.2 

29.9 

934.8 

602.9 

501.7 

1.2 

1479.5 

942.1 

954.0 

-11.8 

1516.4 

1025.7 

985.2 

40.4 

1531.6 

982.1 

998.2 

-16.1 

•H5  GM  TUNGSTEN  SPHEROID  INTO  .32  CM  MILD  STEEL  RT  0 DEGREES  OBLIDUITY 

FIGURE  A- 4 


TABLE  A-5 


• 45  GM  TUNGSTEN  SPHEROID  INTO  .15  CM 
ALUMINUM  AT  45  DEGREES  OBLIQUITY 

CRITICAL  VELOCITY  IS  CONSTRAINED  TO  THE  INTERVAL  FROM  266.4 
TO  410.0  M/S. 


NEI BULL  MODEL 

VC  M/S 

LAMBDA 

BETA 

ERMS 

INITIAL  ESTIMATES 

266.4 

.034665 

.543401 

FINAL  ESTIMATES 

266.4 

.028704 

.623741 

62.9 

STRIKING 

RESIDUAL 

APPROXIMATION 

ERROR 

VELOCITY 

M/S 

VELOCITY 

M/S 

M/S 

M/S 

266.4 

404.8 

.0 

226.5 

187.3 

39.1 

468.8 

297.2 

255.6 

41.6 

480.7 

289.0 

268.2 

20.8 

507.2 

254.2 

296.4 

-42.2 

584.9 

238.4 

379.3 

-140. 9 

907.1 

804.7 

726.9 

77.7 

996.1 

868.4 

823.5 

44.9 

1417.6 

1285.0 

1279.7 

5.4 

1571.2 

1397.8 

1444.7 

-46.9 

?E  A-5 


CRITICAL  VELOCITY  IS  CONSTRAINED  TO  THE  INTERVAL  FROM  482.2 
TO  600.0  M/S. 


ME I BULL 

MODEL 

VC  M/S 

LAMBDA 

BETA 

ERMS 

INITIAL 

ESTIMATES 

482. 2 

.114710 

.390246 

FINAL  ESTIMATES 

581.7 

.344483 

.226746 

43.2 

STRIXING 

RFSIOUAL 

approximation 

ERROR 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

482.2 

.0 

597.7 

291.7 

284.4 

7.3 

666.9 

383.4 

407.4 

-23.9 

695.9 

395.6 

442.0 

-46.4 

880.9 

663.9 

629.7 

34.2 

912.6 

695.9 

659.8 

36.1 

952.5 

727.9 

697.4 

30.5 

1516.4 

1179.0 

1217.7 

-38.7 

1531.6 

1166.2 

1231.7 

-65.5 

1555.1 

1322.8 

1253.3 

69.5 

i 


STRIKING  VELOCITY  M/S 

45  GM  TUNGSTEN  SPHEROID  INTO  .32  CM  RLUMINUM  RT  45  DEGREES  OBLIQUITY 

FIGURE  A-6 


TABLE  A- 7 


.45  GM  TUNGSTEN  SPHEROIO  INTO  .15  CM 
MILO  STEEL  AT  45  OEGREES  OBLIQUITY 

CRITICAL  VELOCITY  IS  CONSTRAINED  TO  THE  INTERVAL  FROM  352.3 
TO  550.0  M/S. 


ME I BULL  MOOEL 

VC  M/S 

LAMBDA 

BETA 

ERM 

INITIAL  ESTIMATES 

352.3 

.273901 

.247218 

FINAL  ESTIMATES 

526.7 

.604598 

.132986 

30. 

STRIKING 

RESIDUAL 

APPROXIMATION 

ERROR 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

352.3 

.0 

541.3 

316.7 

3L3.1 

3.6 

594.1 

408.1 

387.9 

20.2 

613.3 

379.8 

407.9 

-28.2 

623.3 

356.6 

418.0 

-61.3 

658.4 

430.4 

451.4 

-21.0 

662.6 

498.0 

455.3 

42.7 

686.4 

477.9 

477.0 

1.0 

705.6 

501.4 

494.2 

7.2 

718.7 

518.5 

505.8 

12.7 

734.6 

542.2 

519.7 

22.5 

912.6 

685.8 

672.4 

13.4 

923.5 

716.6 

681.7 

34  * 9 

923.5 

681.8 

681.7 

.2 

949.5 

723.3 

703.6 

19.7 

1015.9 

697.1 

759.7 

-62.6 

1508.8 

1178.1 

1175.9 

2.1 

1516.4 

1153.1 

1182.4 

-29.3 

1524.0 

1239.0 

1188.8 

50.2 

1539.5 

1176.2 

1202.0 

-25.8 

SI 


TUNGSTEN  SPHEROID  INTO  -15  CM  MILD  STEEL  ftT  45  DEGREES  OBLIQUITY 

FIGURE  A-7 


TABLE  A- 8 


*45  GM  TUNGSTEN  SPHEROIO  INTO  .32  CM 
MILO  STEEL  AT  45  OEGREES  OBLIQUITV 

CRITICAL  VELOCITY  IS  CONSTRAINED  TO  THE  INTERVAL  FROM  800.1 
TO  890.0  M/S. 


VtE  I BULL  MODEL 

VC  M/S 

LAMBDA 

BETA 

ERMS 

INITIAL  ESTIMATES 

800.1 

.207160 

.200348 

FINAL  ESTIMATES 

877.6 

.396242 

.101630 

58.9 

STRIKING 

RESIDUAL 

APPROXIMATION 

ERROR 

VELOCITY 

VELOCITY 

M/S 

M/S 

M/S 

M/S 

800*1 

.0 

880.9 

317.6 

317.8 

-.2 

955.5 

454.2 

439.9 

14.2 

979.9 

447.4 

460.2 

-12.8 

1516.4 

843.7 

810.0 

33.7 

1516.4 

695.6 

810.0 

-114.5 

1539.5 

903.4 

824.4 

79.0 

STRIKING  VELOCITY  M/S 

GH  TUNGSTEN  SPHEROID  INTO  -32  CM  MILD  STEEL  RT  H5  DEGREES  OBLIQU 

FIGURE  A-8 
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